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Abstract 

Recently, Vatan and Williams utilize a matrix decomposition of SU (2") introduced by Khaneja and Glaser to 
produce CNOT-efficient circuits for arbitrary three-qubit unitary evolutions. In this note, we place the Khaneja 
Glaser Decomposition (KGD) in context as a SU{2") = KAK decomposition by proving that its Cartan invo- 
lution is type AIII, given n > 3. The standard type AIII involution produces the Cosine-Sine Decomposition 
(CSD), a well-known decomposition in numerical linear algebra which may be computed using mature, stable 
algorithms. In the course of our proof that the new decomposition is type AIII, we further establish the follow- 
ing. Khaneja and Glaser allow for a particular degree of freedom, namely the choice of a commutative algebra 
o, in their construction. Let x" be a SWAP gate applied on qubits 1, n. Then /"vx" = ki ak2 is a. KGD for 
a = spanK{x'i'(li)(^-;- 1| " lA'-y- l)(;l)x'i} if and only if v = (x'i*iX'i)(XiO(X?*:2X'[) is aCSD. 

Any fixed-time closed-system evolution of the n-qubit state space may be modelled mathematically by mul- 
tiplication of a state vector |\|/) by some N x N unitary matrix v, where N = 2" throughout. By choice of global 
phase, we may multiply v by det(v)'/^, so that without loss of generality v G SU{N), the Lie group [Kna98, HOI] 
of determinant one unitary matrices. Matrix decompositions are algorithms for factoring matrices. In the context 
of qubit dynamics, such a decomposition would split an evolution v into component subevolutions. 

A matrix decomposition may be proven by explicitly specifying an algorithm that computes it. Alternately, 
several theorems in Lie theory posit factorizations of a group across certain subgroups. These may in essence be 
viewed as meta-decomposition theorems; they often allow many degrees of freedom in the choice of subgroups 
and the group being factored. Without describing the appropriate hypotheses, we mention examples such as the 
global Cartan decomposition G =- exp(p)/r, the Iwasawa decomposition G ~ NAK, and its generalization the 
Langlands decomposition G = UMAK. For G = Gl{n,C) the Lie group of all invertible complex matrices, well- 
known algorithms often exist for the outputs of these theorems. For example, the global Cartan decomposition 
outputs the usual polar decomposition which writes a matrix as a product of a Hermitian and unitary matrix, while 
the Iwasawa decomposition reduces to the QR decomposition. 

The G = KAK metadecomposition theorem has seen several overt and hidden applications in quantum comput- 
ing. A sample output is 5f/(2') = {/;,(e)}{/;;(a)}{/?y(e)},where7?y(e) = cosf|0)(0|+sinf|0)(l|-sinf |0)(1| + 
cos and/?j(a) = e"'"/^|0) (0| +e'"/^|l)(l|. This is the factorization of any Bloch-sphere rotation into rota- 

tions about orthogonal axes. For compact groups such as SU{N), the definitive statement of the G = KAK theorem 
is found in Helgason [HO I , thm8.6,§VII.8]. 

The Khaneja Glaser decomposition (KGD) is constructed by an explicit invocation of the theorem, with G = 
SU (2"). In fact, one formulation [ KGO I , Cor.3] requires two applications of the theorem. In this note, we discuss 
the statement of Corollary 2 ibid. Also, we note that the Cosine Sine decomposition (CSD) [GvL96, pg.77] 
[PW94] of numerical linear algebra is the output of the theorem for one standard choice of inputs for G = SU{2"). 
Per the statement of the abstract, these two matrix decompositions are closely related. Indeed, one results from 
the other after swapping the labels on the first and last qubit. We derive this result in context, using the G ~ KAK 
language. As such, this note is another instance of a Lie theoretic decomposition specializing to matrix analysis. 

The theorem has three inputs, each dependent on the last. The first is the Lie group G with Lie algebra g. 
Should G C Gl{n,C) as a closed subgroup, i.e. G is linear, then g is the set (in fact vector space) of matrix 
logarithms of elements of G. The algebra operation in this case is given by [X,Y] = XY — YX for X,Y G 2- Ths 
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second input is a Cartan involution on g. In the case that G is compact, we take this to mean an R-linear map 
9 : g ^ g is a Lie algebra homomoprhism (9[X,y] = [9X,9y]) and an involution (9^ = ^g-) It is now typical to 
write g p ® t, with p the — 1 eigenspace and t the subalgebra given by the + 1 eigenspace. The final input to the 
theorem is a maximal commutative subalgebra a of p. We now state the G = KAK theorem for compact groups. 
This statement is a consequence of the citation and follows by ignoring the discussion of Weyl group actions ibid. 

Theorem: (Cf. [HO?, thm8.6,§VII.8]) Let U he a connected, compact Lie group with semisimple Lie algebra u, 
let 9 : u ^ u a Cartan involution of u, and let a C p a commutative subalgebra which is maximal with this property. 
For the Lie group exponential, label K = exp(t), A = exp(a). Then 

U=KAK = {kiak2 ; kj e KJ = l,2,a e A } (1) 

Note that the factorization of any given m e t/ is not unique. 

As a brief motivation, we mention how the Singular Value Decomposition (SVD) [GvL96, pg.70] arises as an 
example of the above construction for noncompact G [Kna98, pg.397]. If G = Gl{n,C), then the set of matrix 
logarithms is g = C"^". Take Q{X) = ~X' , so that p the —1 eigenspace coincides with the set of Hermitian 
matrices and i the +1 eigenspace is the set of antiHermitian matrices. Hence K = exp(t) is the unitary group, and 
the Cartan involution has formalized a polar decomposition at the level of g. Now a suitable choice of a would be 
the diagonal matrices, so that Gl{n,C) = KAK = U{n)AU{n) in this case would be the SVD, up to ordering and 
positivity of the diagonal factor In fact, this degree of freedom is accounted for by the Weyl group action in the 
cited theorem. We do not exploit it here. 

The Cosine-Sine Decomposition (CSD) [GvL96, pg.77] may also be viewed as an example of this theorem. 
All possible Cartan involutions 9 of semisimple real Lie algebras, both compact and noncompact, are classified 
[HOI, p.518]. Such Cartan involutions, very loosely analogous to a polar decompositions at the Lie algebra level, 
are less rigid for compact than noncompact algebras. For on a noncompact algebra (e.g. sI(«,M),) all subalgebras 
fixed by Cartan involutions are Lie algebra isomorphic (e.g. to so(n)). In the compact case, this is false. Three 
types of t algebras may arise for Cartan involutions of su(A^). The types are AI, AH, and AIII respectively 
[HOI, pg.518] corresponding to t = so{N), sp{N/2), ands[u{p)®u{q)],p + q=N. The CSD is a SU{N)= KAK 
decomposition whose Cartan involution is of type AIII. 

_ f In/2 



Indeed, take /^/2.«/2 = ( q j 1 = (o'^) i8i/Ar/2- Then the standard Cartan involution in the compact 

V N/2 / 

case of type AIII is 9aiii(^) = In/2.n/2^In/2,n/2 [HOI, pg. 452]. Note that the +1 eigenspace 6aiii is 

«Aiii=s[u(A^/2)®u(A^/2)] = J f ^ ^) ■,X,Y eu{N)MX) + triY)=o\ (2) 



^ y ^ 

There is moreover listed a standard choice of maximal commutative subalgebra. 




aAiii=< J. „ ; r = diag(fo,fi,--- ,fw/2-i) } (3) 



Then the assertion that SU (2") = ^Aiii^Aiii^Ain is the CSD, where ^Taih = exp(tAiii) and Aaiii = exp(aAiii). For 
upon exponentiating, the metadecomposition theorem for these inputs states that any v E SU (N) may be written 



Ml \ / cos(r) -sin(r) \ f U3 
U2 J \ sin(r) cos(r) J \ U4 



(4) 



where uj is unitary for j = 1,2,3,4, the determinant of each factor is one, and T is as in the definition of ttAiii. 
Thus, if y is the namesake cosine-sine matrix of the central factor, we recover the CSD by v = (mi © M2)y("3 © "4)- 
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For the remainder, we use the notations A, K, a, and ? to correspond to the choices of metadecomposition input 
made in the construction of the KGD [KGOl]. Motivated by physical intuition in terms of spin chains, Khaneja 
and Glaser do not expHcitly formulate 9 but rather specify the eigenspace t. This is sufficient, since p and hence 9 
may be recovered from { as a Killing form [HOI, pg.l31] orthogonal complement. The conventions for the KGD 
are then as follows; 

• « = spanR{X®o^y(g)/2,//iv/2«)0^ ■,X,Y <E 5u(A^/2)}, per [KGOl, thm.3] 

• any suitable [KGO 1 , Cor2] maximal commutative a C p 
For the second item, we take %" to be SWAP on qubits 1, «. Label 

a = ^^2.n^{x\{\j){N-j-\\-\N-j-\){j\)x'{\0<j<N/2-\] (5) 

This a is not related to the specific choice of commutative subalgebra of the typical KGD [KGOl, Notations]. 
Rather, it recovers the CSD. 
Proposition: Suppose « > 3 throughout. 

1. Uiii = Xi'6x" 

2. S[U{N /2)®U[N l2)] = t{Kx\ 

3. OAiii = XiaXi 

4. A^^r^Atl 

Thus for V e SU [N], we have v = kyakj a Khaneja Glaser decomposition if and only if X"kjX" — uij-i © U2j € 
S[UiN/2) ®U{N/2)], j = 1 , 2 and x'/a^^ = Y e Aaiii, i.e. if and only if 

XiV'Xi = {rxkitMatiWitx) = («i®«2)y(«3®«4) (6) 
is a Cosine Sine decomposition. 

Proof: For Items 1 and 3, use the standard properties of a SWAP. For the remaining two items, recall that SU [N] 
is a linear algebraic group. Hence the Lie exponential is given by the usual power series for a matrix exponential, 
and generally for any v e SU [N) and X G su(A^) one has exp(vXy^) = v exp(X) y^. Now consider X in each of the 
two Lie algebras and x" = (x")^- ^ 
Although the commutative subalgebra a in the argument above is new, this is not an essential obstacle to 
recovering the KGD from a numerical CSD. For the theory of Weyl group actions demands that all a satisfying 
the hypothesis of the G = KAK theorem are conjugate under some k G K [KGOl, HOI]. Thus, let a ~ X'i^{n)x" be 
the Khaneja Glaser commutative subalgera f)(«) [KGOl, Notations] up to qubit SWAP. Note that specifically 

a = span^l ia^ ® {(f)''^ <S) {(J')''^ ® • • • {(f)''"-^ ® (a-'')®^ ; j ^ Q^x,y,zM e {0, 1},2 < /t < « - 3 } (7) 

Then the theory demands some k € S\U [N /2) © U{N/2)] block diagonal so that ^a^^ ~ flAiii- Indeed, an im- 
portant technique in computing matrix decompositions for two-qubit logic-circuit synthesis is transforming the 
commutative Lie algebra speLn^{i{a'')'^^ Jia^f'^ ,i{a-f'^} into a diagonal Lie algebra [ZVSW03, BM03]. This 
may be accomplished for example by 



1 
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(I /■ 
10/ 
0-1 / 

\ 1 -/ y 



(8) 



Then directly E"^ (a"" » CT^)^ = (g) a^ E"^ (a^ (g) &')E = -o' /2, and £ ' (a^ (g) a'')E = /2 g) a^. Moreover Ha'H = 
H&H'f ^G-\foTH = ^l],^^Q{-iy''\k){j\ the Hadamardmap. Also, if 5 = |0) (0| + !| 1)(1 1, then 5(!0^")5^ = !0>'. 

Thus suppose we take A; = 5 g) g) E^. Then the matrix k may be used to switch commutative subalgebras: 
k~ak^ = {kx'lMn)ikx1)' = span^l (/o>') ® {a')'^ (g • • • (g ; biM, . . . = 0, 1 } = aAiii (9) 
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Thus this block-unitary k S S[U{N/2) ® U{N/2)] allows for translation between the CSD and that KGD where 
the commutative subalgebra is chosen as \){n). 

The CSD has been used for the design of quantum logic circuits directly in [T99, MVBS04, SBM04]. We 
refer to the latter for advice on obtaining a numerical implementation of the CSD. Thus, this document serves to 
translate earlier circuit constructions [V\V()4] in the KGD into more recent works. Similar translations might be 
possible from other applications of the KGD, e.g. in control theory [KGOl]. 
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